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ABSTRACT 

The supercritical disk accretion flow with radiatively driven outflows is studied based on two- 
dimensional radiation- hydrodynamic simulations for a wide range of the mass input rate, M; npu t, 
which is the mass supplied from the outer region to the disk per unit time. The a-prescription is 
adopted for the viscosity. We employ a = 0.5, as well as a = 0.1 for Minput > 3 x 10 2 Le/c 2 and 
only a = 0.5 for M; npu t < 10 2 Le/c 2 , where Le is the Eddington luminosity and c is the speed of 
light. The quasi-steady disk and radiately driven outflows form in the case in which the mass input 
rate highly exceeds the critical rate, M; npu t > 3 x 10 2 Le/c 2 . Then, the disk luminosity as well as 
the kinetic energy output rate by the outflow exceeds the Eddington luminosity. The moderately 
supercritical disk, Mi nput ~ 10 — 10 2 Le/c 2 , exhibits limit-cycle oscillations. The disk luminosity 
goes up and down across the Eddington luminosity, and the radiatively driven outflows intermittently 
appear. The time averaged mass, momentum, and kinetic energy output rates by the outflow as well 
as the disk luminosity increase with increase of the mass input rate, oc M^^'° for a — 0.5 and 
rx M® n p~®' 6 for a = 0.1. Our numerical simulations show that the radiatively driven outflow model for 
the correlation between black hole mass and bulge velocity dispersion proposed by ISilk fc Reed and 
King is successful if Mi nput c 2 /LE ~ a few 10 (a = 0.5) or > a few (a — 0.1). 

black hole physics — galaxies: nuclei — hydrodynamics 



Subject headings: accretion: accretion disks 
— radiative transfer 



1. INTRODUCTION 

The evolution of supermassive black holes (BHs) and 
their host galaxies is a major topic of current inter- 
est. The observations of high-redshift quasars showed 
that t he supermassive BHs had formed in the early uni- 
verse ([Becker et alj l2001t iFan et al.1120011 IShields et all 
l2006h . It has been reported that the gas accretion is 
a sig nificant process for ev olution of the supermassive 
BHs (|Yu fc Tremainej [20021 . If the supercritical accre- 
tion (the mass accretion rate exceeds the critical rate, 
Le/c 2 , where Le is the Eddington luminosity and c is the 
speed of light) is possible, a growth timescale of BHs can 
be much shorter than the Eddington timescale. Thus, 
the supercritical accretion might resolve the problem of 
the formation of supermassive BHs. 

The supercritical accretion onto the central BH is 
also thought to play important roles for the evolu- 
tion of their host galaxies. It has been suggested 
that the correlation between the velocity dispersion of 
the bu lge stars (<Q and th e BH mass (Mbh — re- 
lation: iGebhardt et all 120001 : iFerrarese fc Merrittl 120001 : 
iGreene fc Hoi [2006) , which implies that there is some 
physical link between the evolution of the BHs and the 
host galaxies, is established by the feedback from the 
s upercr i tical a ccretion flows. 

iKingj (|2003f ) suggested that the strong outflow from 
the supercritical accretion flow regulates the evolution of 
the supermassive BH and i ts host galaxy, lea ding to the 
Mbh — c* relation (see also ISilk fc Ree^ll998h . However, 
the momentum and kinetic energy output rates by the 
outflow were assumed without being carefully treated, 



even though they are significant physical quantities in 
t his scenario . 

lUmemural (|2001[ ) suggested that the Mbh — 0* relation 
is built up by mass accretion onto the galactic center 
via t he radiation drag (see also IKawakatu fc Umemural 
|2002[ ) . Since mass accretion is caused by the stellar radi- 
ation of the bulge stars, the BH mass is correlated with 
the bulge mass. Though the mass accretion rate exceeds 
the critical rate in this mechanism, the feedback from the 
supercritical accretion is not taken into consideration. 

Althoug h the standard disk model , which was pro- 
posed by iShakura fc Sunvaevl (|1973f ) . elucidates the 
fundamental properties of luminous compact objects, it 
breaks down in the supercritical accretion regime. In 
this case, a large part of the photons generated inside 
the disk via the viscous process is advected inward and 
swallowed by the BH with accreting matter. This is the 
so-called phot on trap pi ng in the disk accretion flows 
(lOhsuga etall 120021 lOhsuga Mineshige. fc Wataral 
120031 IShimura fc Manmotdl2003h . The radiation pres- 
sure becomes dominant over the gas pressure and drives 
the radiatively driven outflows. The circular motion 
as well as the convection occurs in the disk. They are 
basically multi-dimensional effects. T hese effects are not 
corre ctly treated in the slim disk (|Abramowicz et ail 
1988). since it is a radially one-dimensional model. 

The two-dimensional radiation hydrodynamic 
(RHD) simulations of the supercritical disk ac- 
cretion flows around t he BH s were initiated by 
Eggum. Coroniti. fc Kat d (119881) a nd improved by 
Okudal ([2002D . lOhsuga et all ([2005) for the first time 
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succeeded in reproducing the quasi-ste ady struc t ure o f 
the supercritical disk accretion flows. lOhsugal (|2006l ) 
performed the numerical simulations of the unstable 
disks with the moderately supercritical accretion rate. 
However, the outflow from the supercritical accretion 
flow is poorly understood. Whereas the radiatively 
driven outflow h as extensively been studied by many 
researches (e.g., iBisnovatvi-Kogan fe Blinnikovl 119771 : 
iWatarai fc Fukuej I1999D . they treat the disk as an 
external radiation source and do not solve the radiation 
transfer. 

In this paper, by performing two-dimensional RHD 
simulations of stable and unstable disks, we report the 
feedback (mass, momentum, and energy output rates) 
from the supercritical disk accretion flows around BHs 
for a wide range of the mass input rate, which is the 
mass supplied from the outer region to the disk per unit 
time. 

2. BASIC EQUATIONS AND NUMERICAL METHOD 

Basic equations a nd our numeri c al me thods are de- 
scribe d in detail in lOhsuga et al.1 (|2005l ) and lOhsugal 
(2006). The set of RHD equations including the viscos- 
ity term is solved by an explicit-implicit finite difference 
scheme on the Eulerian grids. Here, we use spherical po- 
lar coordinates (r, 9, ip), where r is the radial distance, 9 
is the polar angle, and <p is the azimuthal angle. The ori- 
gin is set at a central BH. Since we assume axisymmetry 
(i.e., d/dtp = 0), as well as reflection symmetry relative 
to the equatorial plane (with 9 = 7r/2), the computa- 
tional domain can be restricted to one quadrant of the 
meridional plane. The domain consists of spherical shells 
of 3rs < r < 500rs and < 9 < tt/2, and is divided into 
96 x 96 grid cells, where rs is the Schwarzschild radius. 
We describe the gravitational field of the BH in terms 
of pseudo-Newtonian hydrodynamics, in which the grav- 
itational pote ntial is given by -GMbh/^ — ^s) as was 
introduced by iPaczvnskv fc Wiital (|1980n . The flow is 
assumed to be non-self-gravitating. The basic equations 
are the continuity equation, 



g + v.(H = o, 

the equations of motion, 
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the energy equation of the gas, 
de 



dt 



V • (et?) = -pV • v - A-kkB + ckEq + <J> vis , (6) 



and the energy equation of the radiation, 
dE 



where p is the mass density, v = (v ri ve,v v ) is the veloc- 
ity, p is the gas pressure, e is the internal energy density 
of the gas, B is the blackbody intensity, Eo is the ra- 
diation energy density, Fq is the radiative flux, Po is 
the radiation pressure tensor, k is the absorption opac- 
ity, q — (q r , qe , q<p) is the viscous force, and <J> v is is the 
viscous dissipative function. 
The radiation force / rad = (f r , fe) is given by 



/rad — —F 0) 

c 



(8) 



where x{= k+P&t /Trip) is the total opacity with or being 
the Thomson scattering cross-section and m p being the 
proton mass. For the absorption opacity, we consider 
the free-free absorpt ion and the bound-free absorption 
for solar metallicity (iHavash i. HoshL & Sugimotol 119621 : 
iRvbicki fc Lightmanl Il979t ) . In the present study, the 
outflow is mainly accelerated by the radiation force on 
the Thomson scattering, although the radiation force on 
the free-free and bound-free absorption is also taken into 
consideration. 
As the equation of state, we use 



P= (7- !K 



(9) 



where 7 is the specific heat ratio. The temperature of 
the gas, T, can then be calculated from 



P 



pkT 
prrip' 



(10) 



at 



-V-(Eqv) = -V-Fo-Vv : Po+AttkB-ckEq, (7) 



where /cb is the Boltzmann constant and p is the mean 
molecular weight. 

The radiative flux and the radiation stress tensor 
are solv ed under the flux limited dif fusion approxi- 
mation (|Levermore fc Pomraninj Il981f ). so that they 
are expressed in term s of the radiation energy density 
(|Turner fc Stond l200lh . It gives correct relations in the 
optically thick diffusion limit and optically thin stream- 
ing limit, respectively. 

In this paper, we assume that only the r</j-component 
of the viscous stress tensor is nonzero, leading q r — 
qe = 0. This component plays important roles for the 
transport of the angular momentum and heating of the 
disk plasma. In our form, it is proportional to aptotai 
in the optically thick limit and ap in the optically thin 
limit, where a and ptotal are the viscosity parameter and 
the total pressure. Our viscosity model is basically the 
same as the a presc r iption of the viscosity proposed by 
IShakura fc Sunvaevl (j 19731) . 

We start the calculations with a hot, rarefied, and op- 
tically thin atmosphere. There is no cold dense disk ini- 
tially. We assume that matter is continuously injected 
into the computational domain through the outer disk 
boundary (r = 500r s , 0.457T < 9 < 0.5n). The injected 
matter is supposed to have a specific angular momen- 
tum corresponding to the Keplerian angular momentum 
at r = 100rs. We set the injected mass accretion rate 
(mass input rate, -Mmput) so as to be constant at the 
boundary. 

Throughout the present study, we assume Mbh = 
IOM0, 7 = 5/3, and p = 0.5. The a parameter is set to 
be 0.5 and 0.1. The time step is ~ 3 x 10~ 6 s, which is 
determined by min [Ar/c,rA9/c], where Ar and A9 are 
the cell sizes in the radial and polar directions, respec- 
tively. 
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3. RESULTS 
3.1. Highly supercritical case 
In Figure [1] we plot the time evolution of the mass 
accretion rate onto the BH Mbh (black line in the 
top panel), the luminosity -£/ ra d (black line in the bot- 
tom panel), and the mass and kinetic energy output 

rates Af out (red line) and Lkin (blue line) for M; npu t = 
10 3 Le/c 2 , where the viscosity parameter is set to be 0.5. 
The luminosity is evaluated by integrating the radiative 
flux at the outer boundary. It almost equals the luminos- 
ity of the disk. The mass and kinetic energy output rates 
indicate the mass and kinetic energy ejected through the 
outer boundary per unit time by the outflow of v r > v csc , 
where v eS c is the escape velocity. 

Since the matter injected from the outer disk bound- 
ary has angular momentum, it accumulates in the com- 
putational domain without accreting onto the BH in the 
early evolutional phase (t<3 s). The matter starts to 

accrete onto the BH at t ~ 3 s. This critical time 
roughly c oincides with the vis cous timescale (see equa- 
tion [2] in lOhsuga et al.l |2005| ). All the quantities stay 
nearly constant when t>3 s. The total mass contained 
within the computational domain is also constant in this 
phase. Thus, we conclude that the flow is quasi steady 
in this phase. The radiation pressure-dominated thick 
disk forms, and the radiatively driven outflows appear 
above and below t he disk. The q uasi-s teady structure 
is similar to that in lOhsuga et al.l (|2005l ). although they 
employ a = 0.1. 

As shown in Figure [TJ the mass accretion rate is about 
2 orders of magnitude larger than the critical rate, and 
the luminosity stays around 4.8Le in the quasi-steady 
state. We also find that the strong and quasi-steady out- 
flow is generated. The mass output rate is comparable to 
the mass accretion rate. The kinetic energy output rate 
exceeds the Eddington luminosity, whereas it is slightly 
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Fig. 1. — Time evolution of the mass accretion rate onto the 
BH (black line in the top panel), the mass output rate (red line), 
the luminosity (black line in the bottom panel), and kinetic en- 
ergy output rate (blue line) for Afi nput = 10 3 Le/c 2 - The viscosity 
parameter a is set to be 0.5. 



smaller than the luminosity. The energy conversion effi- 
ciency, (L ra d + ikin) /Mbhc 2 , is much smaller than the 
prediction of the standard disk, ~ 0.1. This is because 
a large amount of photons generated inside the disk are 
swallowed by the BH without being radiated away by the 
photon trapping. 

3.2. Moderately supercritical case 

Next, we show the time evolution in the case in which 
the mass input rate moderately exceeds the critical rate 
in Figure O Here, we set the mass input rate to be 
10 2 Le/c 2 and the viscosity parameter to be 0.5. In this 
case, the mass accretion rate drastically varies, although 
the matter is continuously injected into the computa- 
tional domain at constant rate. The mass accretion rate 
suddenly rises from ~ 0.3L E /c 2 to - 300L E /c 2 , and it 
decays after 10 s. Such time variation of the mass ac- 
cretion rate occurs at intervals of ~ 60 s and triggers 
off the luminosity oscillations. The luminosity is around 
0.03Le in the low-luminosity state and 3L E in the high- 
luminosity state. The photon trapping reduces the lumi- 
nosity in the high state. 

Such disk oscillations are induced by the disk 
instability in the radiation pressure-dominated re- 
gion. The disk theory predicts that the radia- 
tion pressure-dominated disk is unstable if the radia- 
tive cooling is dominant o v er the advective cooling 
(iLightman. fe Eardlevl fl97l IShibazaki fe HoshJ H975T 
iKato. Fukue. fe Mineshigd Il998f h Since the advective 
cooling becomes the main cooling process when the mass 
accretion rate is much higher than the critical rate, the 
disk with a highly supercritical accretion rate is sta- 
bilised. In contrast, the moderately supercritical disk 
exhibits limit-cycle oscillations. 

In Figure El we also find that the strong outflow ap- 
pears only in the high-luminosity state. In this state, the 
mass output rate is about 2 orders of magnitude larger 
than the critical rate. The kinetic energy output rate 
exceeds the Eddington luminosity and is comparable to 
the luminosity. In contrast, it is found that L^in is much 
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Fig. 2.— Same as Fig. Q] but for M input = 10 2 L E /c 2 . 
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Fig. 3. — Same as Fig. [T] but for Mj np ut = 10L E /c 2 . 

smaller than L ra d in the low-luminosity state. It im- 
plies that the matter is not effectively accelerated by the 
radiation force when the disk shines at sub-Eddington 
luminosity The radiatively driven outflows are intermit- 
tently produced in the case in which the mass input rate 
moderately exceeds the critical rate. 

The mass ejected from the computational domain per 
unit time (sum of the mass accretion rate onto the BH 
and the mass ejection rate through the outer boundary) 
almost equals the mass input rate on average. Thus, 
simulated limit-cycle behaviour is not a transient phe- 
nomenon. It would continue as long as the mass input 
rate does not change so much. Here, we note that such a 
balance between mass input and ejection rates is for the 
first time reproduced in the present simulatio ns, since 
our in tegration time is much longer than that in lOhsugal 
(2006) and the relatively large viscosity parameter is em- 
ployed. 

Our simulations show that the disk is also unstable 
and exhibits the limit-cycle oscillations in the case of 
-Minput = 10Le/c 2 (see Fig. [3]). As compared with the 
case of Minput = 10 2 £e/c 2 , the bursting behaviour is not 
remarkable. The luminosity, as well as the kinetic energy 
output rate, is comparable to or less than the Eddington 
luminosity even in the high-luminosity state. 

3.3. Feedback 

Based on the results of our numerical simulations for 
a = 0.5, we represent in Figure [4] the time averaged 
mass accretion rate onto the BH (triangles), luminos- 
ity (squares), mass output rate (red circles), kinetic en- 
ergy output rate (blue circles), and momentum output 
rate M out u (green circles) as functions of the mass in- 
put rate. Here, the momentum output rate means the 
radial component of the momentum carried out of the 
computational domain per unit time through the outer 
boundary by the outflow of v r > v csc . The disks with 
Mi nput > 3 x 10 2 Le/c 2 are stable. In contrast, the 
disks exhibit the limit-cycle behaviour in the case of 




Fig. 4. — Time-averaged mass output rate (red circles), the ki- 
netic energy output rate (blue circles), and the momentum output 
rate (green circle) by the outflow for a = 0.5. The mass accretion 
rate onto the BH and the luminosity are plotted with triangles 
and squares. Red, blu e, gree n, and black lines represent the fitting 
formulae of equation (11 ll 14ft . 

We find that all quantities monotonously rise with an 
increase of Mi nput . We present simple fitting formulae 
with the least squares method for our numerical results; 
feedback from the supercritical disk accretion flows in 
the case of a = 0.5 are given by 




0.42, 



2.2, 



1.3, 



1.5. 



(11) 



(12) 



(13) 



(14) 



These fitting formulae are also plotted by red, blue, 
green, and black lines in Figure [4J respectively. Here, 
it is noted that they are valid only in the supercritical 
accretion regime. If the mass input rate is comparable 
to or less than the critical rate, the standard disk, which 
is not accompanied by the outflow, would form. 

Whereas Figured] shows the feedback effects in the case 
of a = 0.5 as functions of the mass input rate, they are 
sensitive to the viscosity parameter. In Figure[5]we show 
the resulting feedback effects for a = 0.1. Here, there are 
no results for M input < 10 2 Le/c 2 . We need to consume 
much time in the numerical simulations for a small mass 
input rate because of the long viscous timescale. In this 
figure, thick solid lines represent the fitting formulae for 
a = 0.1. They are 
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+ 0.62 (15) 
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Fig. 5. — Same as Fig. [4] but for a = 0.1. Thick solid lines rep- 
resent the fitting formulae of equation II15I18H ■ They are extended 
onto the small mass input regime (dotted lines). The fitting for- 
mulae for a = 0.5 are plotted by thin lines. 
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-0.52. (18) 



Though they are derived based on the results for 
A^input > 3 x 10 2 Le/c 2 , we extend them onto the small 
mass input regime (dotted lines). The fitting formulae 
for a — 0.5 are also plotted by thin lines. As shown in 
this figure, we find that feedback effects are more sensi- 
tive to the mass input rate for a = 0.5 than for a = 0.1. 
It is also found that the fitting formulae for the two cases 
intersec t around Mp aut = 3 x 10 2 Le/c 2 . 

Since iKind (|2003l ) proposed that the Mbh — 0* rela- 
tion is established via radiatively driven outflow from the 
supercritical accretion flow onto the supermassive BH in 
the galactic center. The momentum output rate M out K 
is assumed to be ~ Le/c in their model. Thus, our re- 
sults imply that their outflow model is successful as long 
as the normalized mass input rate, Mi nput c 2 / Xe, is a few 
tens in the case of a — 0.5 (see the green line in Fig. 

In the case of a — 0.1, the fitting formula for the 
momentum output rate predicts that the King model is 
successful for Mi nput c 2 / Le ~ a few (see the green dotted 
line in Fig. However, we should note that the fitting 
formulae for a — 0.1 are derived based on the results for 
high mass input rate and there are no numerical results 
for small mass input rate. The fitting formulae might 
overestimate the power of the radiatively driven outflow 
around Mj npu t ~ Le/c 2 . This is because the radiation 
force, which is comparable to the gravity, cannot effi- 
ciently accelerate the outflow. In this case, the normal- 



ized mass input rate for establishment of the Mbh — o~* 
relation would shift towards the higher Mi nput side. We 
need further studies to understand the a dependence in 
detail. 

4. DISCUSSION 
4.1. Formation of supermassive black holes 

Our simulations show that the highly supercritical 
disk is stable and the moderately supercritical disk ex- 
hibits limit-cycle oscillations. The physical mechanism 
for mass supply from the host galaxy onto the galac- 
tic center is one of the most hotly debated issues. Al- 
though its time evolution is not understood yet, the nor- 
malized mass input rate (oc Mi npu t /A/bh ) would grad- 
ually decrease as the seed BH grows into the super- 
massive BH. Thus, our results imply that the seed BH 
grows via the supercritical disk accretion flow in the 
early phase, and subsequently, the BH accretion disk 
evolves into the sub-critical phase by way of the oscil- 
lation phase. The growth timescale via supercritical ac- 
cretion is Mbh/Mbh = 4.5 x 10 6 (M B HC 2 /10 2 J L E )" 1 yr. 

Here, we note that our results would not change that 
much even if we employ a large BH mass (M B h 3> 
10M©), since the radiation force on the Thomson scatter- 
ing mainly supports the disk and accelerates the outflow 
in our simulations. The gas density decreases with an in- 
crease of the BH mass, p oc Mb H > on the condition that 
the normalized mass input rate is kept constant. Hence, 
the radiation force on the Thomson scattering (oc p) is 
dominant over that on the bound-free and free-free ab- 
s orptions (oc p 2 ) ev en in the case of massive BHs. 



lUmemural ([2001D suggested the model for quasar for- 
mation, by which the mass is supplied onto the galactic 
center via the radiation drag, forming the supermassive 
BH in ~ 10 9 yr. Although the mass supply rate onto the 
galactic center highly exceeds the critical rate in their 
model, they assumed the accretion rate onto the BH to 
be the critical rate. The supercritical accretion onto the 
central BH would lead the emergence of quasars in the 
early evolutional phase (-C 10 9 yr). The radiation and 
radiatively driven outflow would affect the stellar evolu- 
tion of the host galaxy. It is not taken into consideration 
in their model. Especially, if the metallicity of the gas is 
much smaller than the local interstellar value, it is known 
that the ultraviolet radiati on plays important roles on 
the star formation activ ity (Hiaman, Rees. fc LoeblH997t 
ISusa fc Umemurall2006[ ). 

4.2. Outflows 

The velocity of the radiatively driven outflow in our 
simulations is 0.1 — 0.3c. It is roughly consistent with 
the results from the observations of quasars. The obser- 
vations of the blueshifted X-ray absorption lines have 
revealed that quasars eject th e highly ionized ma t- 
ter with a velocity of ~ 0.1c (|Pounds et al.l l2003al lbl; 
iReeves. O'Brien, fc Waxcfll200lh . 

Line opacities are not taken into consideration in our 
simulations. The outflow would be enhanced via the radi- 
ation force on lines, unle ss the matter is highly ionized b? 
strong radiation fields fProga. Stone, fc Kallman 200 
IChelouche fc Netzerll2003HProga fc Kallmanll2004r ). The 
line force is thought to be one of the acceleration mech- 
anisms of the outflows in active galactic nuclei (AGNs). 
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The AGNs likely have powerful outflows even for L rac i < 
Lp. 

iProgal ([l999) studied the disk winds driven by line 
force for L ra d < -Le, based on th e stell ar wind model 
given by ICastor. Abbot, fe Kleml (|1975f ). They con- 
cluded that the mass output rate scales with the lumi- 
nosity as Mout cx (I/ ra d/LE) 1_2 ' 5 - Our results do not 
deviate from their result. We derive the relations of 
M out cx (irad/^E) 1 ' 4 from equations (fTTj) and (fT4")) . as 
well as Mont oc (£ ra d/ ^e) 1 ' 6 from equations (fl~5|) and 
(fT5|) . Note that the disk is treated as the source of radia- 
tion and mass without solving its structure in their work. 
In contrast, both the disk and the outflow are solved in 
our simulations. In addition, the central BH swallows a 
large amount of matter in our work. 

The magnetic effects might also enhance the out- 
flow. It has been reported by magnetohydrody- 
namic (MHD) simulations that the magnetized accre- 
tion disks produce the outflows, although these simu- 
lations focus on the rad i atively inefficient accretion flow 
(lHawlev fc Balbusll2002t iKudoh. Matsumoto. fc Shibatal 
120021 iKato. Mineshige. fc Shibatall2004fj ~ 

4.3. Future work 

In this paper we investigate the feedback from the su- 
percritical disk accretion flows around BHs as functions 
of the mass input rate. Although we employ merely two 
viscosity parameters, a = 0.5 and 0.1, the feedback ef- 
fects might be more sensitive to the viscosity parameter 
than to the mass input rate. In addition, we do not 
present results for Mi npu t < 10 2 Le/c 2 in the case of 
a = 0.1. Further long integration time is required to 
simulate the accretion flows with a small viscosity pa- 
rameter and/or a small mass input rate because of the 
long viscous timescale. Such numerical simulations are 
left as future work. 

More importantly, we need radiation MHD simula- 
tions, since the sour ce of disk viscosity is likely to be 
of magnetic origin (IStone. Pringle. fc Begelmanl 1999 : 



work. 



5. CONCLUSIONS 



Stone fc Prinriel 2001: Hawlcv. Balbus. fc Stone! 2001 



Machida. Matsumoto. fc Mineshigel 20011 ). A more phys 
ical treatment of the viscosity might lead to the difference 
from the a-disks. Whereas local radiatio n MHD simu- 
lation s have been performed recently (e.g.. iTurner et all 
2003), global simulations should be explored in future 



By performing the two-dimensional RHD simulations, 
we study the feedback (mass, momentum, and energy 
output rates) from the supercritical disk accretion flows 
around BHs for the wide range of the mass input rate, 
-Minput = 10 — 3 x 10 3 £e/c 2 . The adopted viscosity pa- 
rameter a is 0.5 and 0.1 for Mi nput > 3 x 10 2 Le/c 2 and 
only 0.5 for M input < 10 2 Le/c 2 . The detailed study for 
a-dependence is left as future work. We summarize our 
results as follows. 

1. In the case of high mass input rate, Afj npu t > 3 x 
10 2 Le/c 2 , the quasi-steady disk forms and the outflows 
driven by the radiation force appear above and below the 
disk. Both energy output rates by the radiation and the 
outflow exceed the Eddington luminosity. 

2. The disk with moderately supercritical rate, 
Minput ^10— 10 2 Le/c 2 , exhibits limit-cycle oscillations, 
leading to intermittent outflows. The strong outflow is 
generated only in the high-luminosity state, in which the 
energy output rates by the radiation and the outflow are 
comparable to or exceed the Eddington luminosity. 

3. The time-averaged mass, momentum, and kinetic 
energy output rates by the outflow, as well as the disk 



»V0.7-1.0 



for 



luminosity, scale with mass input rate as oc M inpu . 
a = 0.5. In the case of a = 0.1, we find that they are 
proportional to M^~° 6 . 

4. Our numerical simulations show that the radiatively 
driven outflow model for the A/bh — o~* relation proposed 
by Silk & Rees and King is successful, as long as the 
normalized mass input rate Mj nput c 2 /LE is around a few 
10 (a = 0.5) or a few (a = 0.1). 
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